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Abstract. The recent proposal by Kellermann (1993) to 
use the angular size/redshift relation of compact radio 
sources to obtain information on the value of cosmolog- 
ical parameters has provoked considerable discussion in 
the astronomical community. I present here a careful re- 
analysis of the angular size/redshift data of radio jets used 
by Kellermann (1993), using a Kolmogorov-Smirnov test 
to check the independence of the distribution of linear jet 
sizes from redshift under the assumption that no evolu- 
tion is present in the sources. The test takes into account 
biasing of the sample due to the limited resolution. 

While the data are compatible with an Einstein- 
de Sitter Universe (P = 0.2), they are not sufficient to 
rule out a significant part of the (Slo, Ao)-plane. 
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1. Introduction 

Recently Kellermann (1993) made the interesting proposal 
to use the angular size/redshift relation of compact radio 
sources to obtain information on the cosmological decel- 
eration parameter qo, which in the case of a vanishing 
cosmological constant equals Slo/2. 

While earlier attempts to use the angular size/redshift 
relation for high redshift optical or radio sources have been 
faced with severe problems due to the difficulty of pre- 
cisely defining and measuring angular diameters, as well 
as excluding evolutionary effects, compact radio sources 
seemingly offer a way around these obstacles. 

Kellermann (1993) concluded that the data are con- 
sistent with an Einstein-de Sitter Universe and that they 
provide direct evidence that the density parameter does 
not differ significantly from unity. However, this conclu- 
sion is based only on visual comparison of the binned data 
with the theoretical angular size/redshift relation for mod- 
els with a vanishing cosmological constant. 



In this letter a careful re-analysis of the original un- 
binned data is presented, without a priori restrictions on 
the cosmological model, using a Kolmogorov-Smirnov test 
and taking into account biasing due to the limited angular 
resolution in the sample. 



2. The sample 

For his analysis Kellermann (1993) used a sample of 82 
compact radio sources taken from the existing literature 
which met his selection criteria discussed below. A publi- 
cation of the data of the complete sample, including ref- 
erences, is in preparation (Kellermann 1994). I will here 
give a short overview of the selection criteria as well as of 
the definition and measurement of the angular size. For a 
complete description see Kellermann (1993). 

In order to confine the sample to a homogeneous group 
with similar physical properties, only sources with a lumi- 
nosity greater than 10 24 W Hz -1 (at A = 6 cm) have 
been used. This is roughly the luminosity which separates 
radio-loud from radio-quiet quasars. Furthermore, Keller- 
mann (1993) excluded all known BL Lac objects, compact 
central components of extended double- lob ed quasars, and 
peaked spectrum sources. 

In addition, I excluded from the sample not only un- 
resolved sources but also sources with redshifts below 
z = 0.03, for which the Doppler redshift due to peculiar 
motion may not be negligible. This reduced the number 
of objects from 82 to 76. 

Kellermann defined the angular size 6 of a source as the 
distance between the core and the most distant component 
whose peak brightness exceeds 2% of the core brightness. 
This definition of the angular size should keep any effect 
of observing at a fixed frequency small. All sizes were de- 
termined from the published contour diagrams based on 
VLBI observations at 6 cm, a resolution 6o of about 1 
milliarcsecond, and a dynamic range of at least 1:100. 
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Fig. 1. Angular size/redshift diagram for the complete sample 
of compact radio sources used in this analysis. 



3. The statistical method 

My aim is, within the framework of standard Friedmann- 
Lemaitre cosmology, to obtain information on the cosmo- 
logical parameters Slo (density parameter) and Ao (nor- 
malised cosmological constant, Ao = A/(3Hq)) from the 
distribution of angular jet sizes 6 in dependence on the 
redshift z. 

Kellermann (1993) collected the data into redshift 
bins, calculated the mean angular size for each bin, and 
than compared visually the binned data to the expected 
angular size/redshift relation for Ao = models. 

While this approach makes sense to get a first overview 
over the data and to demonstrate in an illustrative way 
how the proposed test works, it is clear that a more sound 
statistical test is needed to obtain a 'best fit' world model 
or to reject world models incompatible with the data. In 
addition, no a priori restriction (like setting Ao = 0) should 
be made for the cosmological model. 

Binning the data introduces several problems into the 
analysis (see, e.g., Press et al. 1986). Since each bin rep- 
resents a subsample of the data by just one quantity (the 
mean), information on the inner distribution of the data is 
lost. The mean of the data is in most cases (including this 
one) not representative of the value of the considered func- 
tion at the centre of the bin. Finally, binning introduces 
free parameters (size and positions of the bins) into the 
analysis, on which the results obviously depend in some 
difficult to control way. The latter is especially true if the 
number of bins and the number of objects per bin is small. 



data without binning for the analysis. Fig. 1 clearly shows 
that the distribution of the angular size 6 is not a nor- 
mal (Gaussian) one. There is some concentration towards 
small 6 and a significant tail at high 6. Thus, a standard 
X 2 approach cannot be used to fit a 6(z) relation to the 
data. 

Assuming that the radio jets show no evolution with 
redshift is equivalent to assuming that the distribution of 
linear jet sizes L' is independent of the redshift. It seems 
natural to test this assumption by means of a Kolmogorov- 
Smirnov test. To do this, I split the sample into a low 
redshift (z < z cut ) and a high redshift (z > z cut ) sub- 
sample, and check whether the distribution of the linear 
sizes (for the considered world model) is the same for both 
subsamples. 

The linear size of the jets is obtained from 



L' = D(H ,Q ,X ;z) 



(1) 



where D = rR is the angular size distance (see, e.g., Berry 
1989), which depends on the Hubble parameter Ho, the 
world model (Slo, Ao) and the redshift z of the source. Since 
I am only interested in the distribution of the jet sizes and 
not in their absolute values, I use only dimensionless linear 
sizes defined as 



L' 



Ho 



Dt 



(2) 



which do not depend on Ho since D oc Hq . 

While the sample is (hopefully) complete with respect 
to 6 down to the resolution limit 6o , it is not complete with 
respect to the linear sizes L. Sources with linear sizes 



L < D(tto,\o;z)6o 



(3) 



are not resolved, a limit obviously depending on the red- 
shift and the world model. For each world model the an- 
gular size distance has a certain maximum value D max for 
the considered sample. I denote 



Lo — D max (Qo, Ao) do 



(4) 



as the linear resolution limit for the world model (Slo, Ao). 
It is clear that the sample is incomplete for L < Lo- 
Sources smaller than the linear resolution limit are not 
resolvable at all redshifts. 

For each considered world model I must therefore first 
calculate D max and Lo, and then exclude all objects with 
L < Lo to avoid biasing due to the limited resolution. 

The Kolmogorov-Smirnov test applied to the two sub- 
samples gives the probability P that the subsamples are 
drawn from the same parent distribution. This is always 
an upper limit to the probability that the distribution of 
L is independent of the redshift. 




icant part of the (Slo, Ao) plane at a sufficient significance 
level (P < O.Of). Note that the white area (P = 0) in 
Fig. 2 corresponds to so-called bounce models (Stabell & 
Refsdal f 966) which are rejected by the mere existence of 
redshifts up to z = 4. 

The different result of Kellermann (1993) is mainly an 
artefact of the binning scheme and of the a priori restric- 
tion to Ao = models. Due to the large intrinsic scat- 
ter of the sizes of compact radio sources, their angular 
size/redshift relation does not provide us with any evi- 
dence that the cosmological constant is close to zero or 
that the density parameter is close to unity. 
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Fig. 2. Probability P(Qo,Ao) that the distribution of linear 
jet sizes is independent of the redshift. White corresponds to 
P = 0. The thick white line, the thin white line, the thin 
black line and the thick black line correspond to P = 0.99, 0.9, 
0.5, and 0.1, respectively. The straight black line marks flat 
universes (k = 0). 

4. Results and discussion 
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I have tested 300 x 300 Friedmann-Lemaitre models in the 
range Slo = • • • 10 and Ao = —5 • • • + 5, using the method 
described above. In order to keep the number of objects 
in both subsamples and for all tested world models above 
20, which is necessary to obtain reliable results from the 
Kolmogorov-Smirnov test (Press et al. 1986), the sample 
is split at z cut = 0.75. In Fig. 2 the resulting probability 
distribution is shown. 

The highest probability (P > 0.99) is found for Slo = 
1.0 •• -3.0 and Ao = 1.5 •• -2.5. This should, however, not 
be interpreted as a "best fit" , since high probabilities in 
KS-tests usually occur just by chance, while the true value 
of the KS-test is the rejection of models with low proba- 
bilities (P < 0.01 is a conservative significance level). 

Marginal evidence (P < 0.1) can be found for the ex- 
istence of (non-baryonic) dark matter, since a minimum 
value of Slo = 0.2 is found for P > 0.1, which is larger than 
the current limit for baryonic matter as obtained from 
nucleosynthesis (see, e.g., Krauss 1994). If one assumes 
Ao = 0, a minimum of Slo = 0.53 is required for P > 0.1. 
If one assumes a flat universe (k = 0, i.e. Slo + Ao = 1), 
Sl = 0.67 is required for P > 0.1. 

The data are compatible with the Einstein-de Sitter 
case Slo = 1 and Ao = at a moderate probability of 
P = 0.2. 
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